Microbial GWAS studies revealing combinations of Omicron RBD mutations existed and may contribute to antibody evasion and ACE2 binding (preprint)

Since Omicron variant of SARS-CoV-2 was first detected in South Africa (SA), it has now dominated in United Kingdom (UK) of Europe and United State (USA) of North America. A prominent feature of this variant is the gathering of spike protein mutations, in particularly at the receptor binding domain (RBD). These RBD mutations essentially contribute to antibody resistance of current immune approaches. During global spillover, combinations of RBD mutations may exist and synergistically contribute to antibody resistance in fact. Using three geographic-stratified genome wide association studies (GWAS), we observed that RBD combinations exhibited a geographic pattern and genetical associated, such as five common mutations in both UK and USA Omicron, six or two specific mutations in UK or USA Omicron. Although the UK specific RBD mutations can be further classified into two separated sub-groups of combination based on linkage disequilibrium analysis. Functional analysis indicated that the common RBD combinations (fold change, -11.59) alongside UK or USA specific mutations significantly reduced neutralization (fold change, -38.72, -18.11). As RBD overlaps with angiotensin converting enzyme 2(ACE2) binding motif, protein-protein contact analysis indicated that the common RBD mutations enhanced ACE2 binding accessibility and were further strengthened by UK or USA-specific RBD mutations. Spatiotemporal evolution analysis indicated that UK-specific RBD mutations largely contribute to global spillover. Collectively, we have provided genetic evidence of RBD combinations and estimated their effects on antibody evasion and ACE2 binding accessibility.

Tracing two causative SNPs reveals SARS-CoV-2 transmission in North America population (preprint)

During the COVID-19 pandemic, precisely tracing the route of the SARS-CoV-2 transmission in human population remains challenging. Because this RNA virus can mutate massively without a specifically tracing maker. Herein, using a geographic stratified genome-wide association study (GWAS) of 2599 full-genome sequences, we identified that two SNPs (i.e., 1059.C>T and 25563.G>T) of linkage disequilibrium were presented in approximately half of North America SARS-CoV-2 population (p = 2.44 x 10-212 and p = 2.98 x 10-261), resulting two missense mutations (i.e., Thr 265 Ile and Gln 57 His) in ORF1ab and ORF3a, respectively. Interestingly, these two SNPs exclusively occurred in the North America dominated clade 1, accumulated during mid to late March, 2020. We did not find any of these two SNPs by retrospectively tracing the two SNPs in bat and pangolin related SARS-CoV-2 and human SARS-CoV-2 from the first epicenter Wuhan or other regions of China mainland. This suggested that the SARS-CoV-2 population of Chinese mainland were different from the prevalent strains of North America. Time-dependently, we found that these two SNPs first occurred in Europe SARS-CoV-2 (26-Feb-2020) which was 3 days early than the occurring date of North America isolates and 17 days early for Asia isolates (Taiwan China dominated). Collectively, this population genetic analysis highlights a well-confidential transmission route of the North America isolates and the two SNPs we newly identified are possibly novel diagnosable or druggable targets for surveillance and treatment.